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Abstract.  Recently,  Parida  and  Gupta  [J.  Comp.  Appl.  Math.  206  (2007), 
873-877]  used  Rail’s  recurrence  relations  approach  (from  1961)  to  approximate 
roots  of  nonlinear  equations,  by  developing  several  methods,  the  latest  of  which 
is  free  of  second  derivative  and  it  is  of  third  order.  In  this  paper,  we  use  an 
idea  of  Kou  and  Li  [Appl.  Math.  Comp.  187  (2007),  1027-1032]  and  modify 
the  approach  of  Parida  and  Gupta,  obtaining  yet  another  third-order  method 
to  approximate  a  solution  of  a  nonlinear  equation  in  a  Banach  space.  We  give 
several  applications  to  our  method. 


1.  Introduction 

Newton’s  method  and  its  variants  are  used  to  solve  nonlinear  operator  equations 
F(x)  =  0  or  systems  of  nonlinear  equations.  The  convergence  of  these  methods 
was  established  using  Kantorovich  theorem  (see  e.g.  [1,  2,  3]).  The  convergence 
of  the  sequences  obtained  by  these  methods  in  Banach  spaces  is  derived  from  the 
convergence  of  majorizing  sequences  (see  [4]  and  references  therein).  Rail  [5]  has 
suggested  a  different  approach  for  the  convergence  of  these  methods,  based  on 
recurrence  relations.  Parida  [6],  and  Parida  and  Gupta  [7]  used  this  idea  for  several 
third-order  methods  (see  also  the  work  of  Candela  and  Marquina  [8,  9],  Ezquerro 
and  Hernandez  [10],  and  Gutierrez  and  Hernandez  [11,  12]). 

Here  we  apply  the  idea  to  the  third-order  method  free  of  second  derivative  pro¬ 
posed  by  Kou  and  Li  [13].  They  developed  a  family  of  methods  for  the  solution  of 
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Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


a  nonlinear  equation  f(x)  =  0  as  follows 


Vn 


xn 


n  f{Xn ) 

r(xny 


*^n+ 1 


f{yn)  +  (92  +  9  -  l)f(xn) 

e2f'(xn ) 


(1.1) 


It  turns  out  that  this  method  is  of  third  order  when  approximating  a  simple  root. 


2.  Recurrence  relations 

In  this  section,  we  discuss  a  third-order  method  for  solving  nonlinear  operator 
equations 

F(x)  =  0,  (2.1) 

where  F  :  Q  C  X  — >■  T  is  a  nonlinear  operator  on  an  open  convex  subset  fl  of  a 
Banach  space  X  with  values  in  a  Banach  space  Y .  The  third-order  method  [13]  is 
defined  as  follows: 


Dn  =  xn  -  9F'{xn)  1F(xn), 

xn+i  =  xn-°  +®2 — -F\xn)~1F(xn)  -  ^F'(xn)~1F(yn). 

This  family  uses  two  evaluations  of  F  and  one  evaluation  of  F' .  In  [7]  they  discuss  a 
third-order  method  requiring  one  evaluation  of  F  and  two  evaluations  of  F' .  Several 
choices  of  9  were  suggested  in  [14]  and  [15]. 

Let  F  be  a  twice  Frechet  differentiable  in  f l  and  BL(Y ,  X)  be  the  set  of  bounded 
linear  operators  from  Y  into  X.  Let  us  assume  that  Tq  =  F'(x o)-1  G  BL(Y,  X) 
exists  at  some  xq  el!  and  the  following  conditions  hold: 

(!)  \\F'(x)  -  F'{y)\\  <  ki\\x  —  y\\,  x,y  G  Q, 

(2)  \\F"(x)\\<M,  xen, 

(3)  ||r0||  < 

(4)  |jr0F(a;o)||  <  r). 

Let  us  also  denote 

a  =  ki/3r], 

a  =  - ~q2 - .  (2-d> 

7  =  fPv- 

Now,  we  define  the  sequences 


a0 

Ora+l 

bn+ 1 
9n+ 1 


6o  —  1,  do  —  a  +  J,  6_i  —  0, 

1  aandn 


an+il3riCn, 

\92  +  9  -l\  1 

- Q2 - &"+1  + 


|1  -  6»|C7„  + 


M 


e2b2n+ 1 


(2.4) 


2 


where 


Cn  =  YKn  +  *ll0l  b^Kn  +  yl02  -  Mbl, 

T .  |0  +  1|(0-1)2  +  |1-0|  .  M  pu2 

Kn  =  - - - J2 - -  K  +  —an(3b2nV. 


Note  that  we  can  rewrite  dn+i  also  in  the  form 


dn-\-\  Qtbn- )_i  ”/CLn-\-lbn-\-h  (2.T) 

or,  equivalently,  as 

dn+ 1  =  dobn+i  +  7&„+i  (a„+i&„+i  —  1) . 

The  polynomials  Cn  and  Kn  can  be  rewritten  as 

Cn  =  (Pq  +  PlO-nbn  +  P2dnbV)  bn> 

dP-n  ~  (Qo  "b  QlCLn.bn )  bn. 

Lemma  1.  Under  the  previous  assumptions,  we  prove  the  following: 

(In)  ||r„||  =  ||T1,(a;n)~1||  <  an(3, 

(I In)  \\UnF(xn)\\<bnri, 

(IHn)  ll^n+l  Xn  ||  ^  dnT], 

(IVn)  \\xn+i  -  y„  II  <  (dn  +  2Kn_i  +  9bn)r). 

Proof.  We  use  induction  to  prove  the  above  claims.  Notice  that  (Jo)  and  (/Jo) 
follow  immediately  from  the  assumptions.  To  prove  ( IIIq ),  we  start  with  the  first 
substep  of  (2.2), 

F(yo)  =  F(y0)  -  0F(xo)  -  F'(x0)(yo  -  x0).  (2.8) 

This  can  be  written  as 

F(yo)  =  (l- 9)F(x0)  +  F(y0)  -  F(x0)  -  F'(x0)(y0- x0)  .  . 

=  (1  -  9)F(x 0)  +  fg  F"(x o  +  t(y0  -  x0))(l  -  t)dt(y0  -  x0)2. 

Now  multiply  by  T0  and  use  the  assumptions,  it  follows  that 

„„  . .  M  O  Q  , 


so  that 


||r0F(t/o)||<|l-%+  y^V, 


11*1 -*o||  <  162 +  °  1|||r0F(xo)||  +  ^||r0F(W))|| 


^  / 102  +  0  —  1|  +  |1  —  0|  M  \ 

<  y - 02 - +  y^J  9  =  d°9, 

and  (IIIq)  holds. 

We  have 

*i  -  yo  =  -  ^  +  1  — — T0F(a;o)  -  ^T0F(y0), 

so  that  it  follows  from  (2.10)  that 

ll*i -Stoll  <  |g  +  1|02~1)2|iro^o)ll  +  yliroFMII 

^  /|0  +  1|(0-1)2  +  |1-0|  Ma  \ 

<  [[ - -  gi - - [  +  YP9)9  =  doV 

<  (do  +  2K_i  +  9b0)rj, 


and  (IV o)  also  holds. 

Following  an  inductive  procedure  and  assuming  xn  €  Q  and  aandn  <  1,  if  xn+i  €  0, 
we  have 


\\I  -TnF'(xn+1)\\  <  ||rn||  HF'^a’n)  —  F'(xn+i)\\  <  aandn  <  1.  (2.14) 

Now,  we  note  that 

rn+1  [I  -  (F'(xn)  -  F'(xn+1))  rn]  =  r„.  (2.15) 

Then  Tn+i  is  defined  and 

||r„||  ^  an/3 


l|r„+i||  < 


< 


1  -  ||rn||||^(*n)  -  U,(a:r!+1)||  1  -  aand, 


=  an+ 1/3.  (2.16) 


Hence,  by  induction  (2.16)  holds  for  all  n.  This  proves  condition  (/„). 

Using  the  first  step  of  (2.2),  we  have 

F(yn)  =  F(yn)  -  0F(xn)  -  F'(xn)(yn  -  xn) 

=  (1  -  9)F(xn)  +  F(yn)  -  F(xn)  -  F'(xn)(yn  -  xn)  (2.17) 

=  (1  -  9)F(xn)  +  /q  F”( xn  +  t(yn  -  ®n))(l  -  t)dt(yn  -  xn)2. 

Now  subtract  the  first  step  of  (2.2)  from  the  second  to  get 

F'( xn)(xn+1  -  yn)  =  - - °  q2  d+1F(xn )  -  ^ F(yn ).  (2.18) 

Using  (2.17)  on  the  identity 

F(xn+ 1)  =  F'(xn)(xn+l  -  yn)  +  F(yn)  +  [F'(yn)  -  F'(xn)\(xn+1  -  yn) 

+F(xn+ 1)  -  F(yn)  -  F'(yn)( xn+1  -  yn) 

(2.19) 

and 

F(xn+\)-F(yn)-F'  (yn)(xn+i-yn)  =  /  F"  (yn-\-t(xn+i-yn))(l-t)dt(xn+1-yn)2 , 

Jo 

(2.20) 


we  have 

F(xn+ 1) 


]^p()-F(xn)  -  ^ F(yn ) 


+  f  F"(xn  +  t(yn  -  xn))(l  -  t)dt(yn  -  xn)2 

(2.21) 

+  /  F"(yn  +  t(xn+ 1  -  yn ) ) ( 1  -  t)dt(xn+i  -  yn )2 
Jo 


+  [F'(yn)  -  F'(xn)](xn+ 1  -  2/n), 


or 


q2  _  ^  F 

F(xn+ 1)  =  ^  J  F" (xn  +  t(yn  -  xn))(l  -  t)dt(yn  -  xn)2 

+  f  F" (yn  +  t(xn+ 1  -  yn))(l  -  t)dt(a:n+i  -  yn)2  (2.22) 
Jo 


+  [F'(yn)  ~  F'(xn)](xn+ 1  -  2/n). 
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Hence, 


\\F(Xn+1)\\  <  ^^A\\yn-Xnf  +  ^i\\Xn+1-ynf 


From  (2.17)  we  get 


+h\\yn  -  xn\\\\xn+i  -  yn\\. 


M 


llrnF1^)!!  <  |1  -  0\ bnV+  anp\\yn  -  xn\\2, 


(2.23) 


(2.24) 


so  that  since  \\yn  —  xn\\  =  \\6rnF(xn)\\  <  \0\bnrj,  and  combining  (2.24)  with  (2.18), 
we  have 

ii  n  ^  —  0  +  1|  +  |1  —  9\  M  on  ..2 

ll^n+l  —  Vn  ||  ^  ~q2  ^ nV  +  ~^Q2anP\\yn  ~  xn  || 

<  l',  +  1K<i-()l>a  +  |1-0|  M  +  =  Knl.  (2.25) 


Hence, 


where  Cn  is  given  by  (2.5). 
Therefore 


||F(xn+i)||  <  Cnrj2, 


||rn+1F(a:n+1)||  <  ||rn+1||||F(a;n+1)|| 

^2  Qn+lfiCnT] 

So,  by  induction  (2.27)  holds  for  all  n.  This  proves  condition  (//„). 
Using  (2.17)  with  xn  and  yn  replaced  by  xn+±  and  yn+i  respectively, 

\\F{yn+1)\\  <  \l-6\ Cny2 +A^92b2n+1ij2, 


so  that 

||r„+iF(yn+1)||  <  ||rn+1||||F(yn+1)||  <  an+1  fiy2 

Using  (2.27)  and  (2.29),  we  therefore  have 
\d2  +  6-  1 


\i-e\cn  +  fe2b2n+1 


(2.26) 

(2.27) 

(2.28) 
(2.29) 


\\xn+2  ~  Xn+i\\  < 


92 


T„+iF(arn+i||  +  ^2  l|rn+iF(j/n+i)|| 


^  \e2  +  e-  ii  i  „  2 

<  - - bn+iV  +  -^an+1(3ri 

\o2  +  8-i\i  ,  1 

- ff2 - 6«+l  +  ff2an+lPV 


\i-o\cn  +  fe2b2n+1 

|1  —  $1  Cn  + 


V 

(2.30) 


—  dn+ 1  ^7 , 

Hence,  by  induction,  this  inequality  holds  for  all  n.  This  proves  condition  (///„). 
Since 

{9  +  1)(6  —  l)2 
9 2 


Xn+l  l/n+1  —  yn  ^n+1  T 


-TnF(xn 


1 


~Q2  TnF{yn)  +  9Tn+1F(xn+1), 
we  have  from  (2.24),  (2. 25), (2. 27)  that 

ll^ra+i  -  Un+ 1||  <  (2 Kn  +  8bn+i)rj. 


(2.31) 

(2.32) 


5 


(2.33) 


Hence  we  have 


\\%n+2  -  Vn+lW  <  \\xn+2  -  Xn+i\\  +  \\xn+i  -  yn+i\\  (2.34) 

<  dn+ir]  +  (2Kn  +  6bn+i)r]  (2.35) 

=  (rfn+i  +  2ivn  +  0&„_|_i)r;.  (2.36) 

Hence,  by  induction,  condition  (IVn)  holds  for  all  n.  This  proves  condition  (IVn). 

a 


3.  Convergence  Analysis 


In  this  section,  we  shall  establish  the  convergence  of  our  third-order  method 
(2.2).  To  this  end,  we  have  to  prove  the  convergence  of  the  sequence  xn  defined  in 
a  Banach  space  or,  which  is  the  same,  to  prove  that  dn  is  a  Cauchy  sequence  and 
that  the  following  assumptions  hold: 

(1)  %n  U, 

(2)  aandn  <1,  n  €  N. 

The  next  two  lemmas  will  show  the  Cauchy  property  for  the  sequence  dn. 

Lemma  2.  Assume  that  Xo  is  chosen  so  as  to  satisfy  0  <  do  <  1/a,  that  is, 
a(a  +  7)  <  1,  where  a  and  7  are  given  by  (2.3).  Then,  the  sequence  an  >  0  is 
increasing,  as  n  increases. 


Proof.  We  show  now  that  all  the  involved  sequences  are  positive.  Under  the  imposed 
conditions,  we  see  that  ao,  bo,  do,  Co,  Kq  are  all  positive,  and  also  that  1  —  aaodo  >  0. 
Assume,  now,  that  all  a,  ,  6,,  di,  Ci,  K, ,  and  1  —  acqdj  are  positive,  for  *  =  0, 1, ... ,  n. 

Since  Cn  >  0  and  bn+ 1  =  an+i(3r]Cn,  it  follows  that  an+i,6n+i  have  the  same 
sign,  and  so  an+ibn+i  >  0.  Further,  from  dn+ 1  =  6ra+i(a  +  7an+i&n+i),  we  get 
that  dn+ 1  has  the  same  sign  as  bn+ 1,  and  so,  all  three  terms  an+ 1,  bn+ 1,  dn+ 1  share 
the  same  sign. 

By  absurd,  we  suppose  that  the  implied  sign  is  negative.  Then  dn+dn+ 1  <  dn, 
and  so, 

1  -  aan(dn  +  dn+i)  >  1  -  aandn, 

which  renders 

1  ,  _  1  aan(dn  dn+\) 

1  —  aan+ia„_|_i  — - - - - - >  1, 

1  -  aandn 

which  implies  aan+idn+i  <  0,  but  that  is  impossible  since  an+i,  dn+\  have  the  same 
sign  and  a  >  0. 


Next,  since  an+ 1  = 


1  —  aa„  dr, 


then 


1 


,  _  -  ,  1  1 

dn  —  I 

CL  \  CLn  CLn-\-l 


and  so,  by  telescoping,  we  get 

n—  1 


L '  n 


i= 0 


1 

a  \a0 


,  where  ao  =  1. 


This  will  render 


n— 1 


1  —  a  di 


i= 0 
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n— 1 

Certainly,  since  a  >  0,  di  >  0,  for  all  i  >  0,  then  a  di  increases  as  n  increases, 

i= o 


n—  1 

and  so,  1  —  a  di  decreases  as  n  increases,  which  implies  that  the  reciprocal, 
i— 0 

namely,  an  is  an  increasing  sequence,  and  consequently,  an  >  «o  =  1.  □ 


We  define  the  sequence  cn  =  anbn.  Then  the  sequences  {an},  {&„},  {cn},  {dn} 
can  be  rewritten  as 


«n+ 1 
bn+1 


^n+1 

dn+ i 


1  -  a  (acn  +  7 c2 ) 7 
/3?7&ncra  (P0  +  Pic„  +  P2c2n) 

l-  a  ( acn  +  7c2 ) 
pr]c2n  (Pq  +  PjCn  +  P2c2n) 

[1  -  a(ac„ +7C2)]2 

a&„+i  +  7an+i62+1 

Pv bncn  (P0  +P1cn  +  P2c2n ) 

- - - - - - - - (a  +  7Cn+i)  • 

1  -  a{acn  +  7C„) 


1  -  aandn 
an+1  Pr]Cn  = 


Gn+l^n+l  — 


That  the  sequence  {cn}  is  a  decreasing  sequence  under  the  assumption  that 
ai&i  <  1  can  be  proved  by  using  the  mathematical  induction.  It  is  obvious  that 
Ci  =  aibi  <  1  =  Cq.  Assuming  that  cn  <  cn_i  for  some  n  >  0,  we  have 


_  Pv cl  (P0  +  Picn  +  P2c2) 

[1  —  a  (acn  +  7cn)]2 
.  Pvcl-l  (Po  +  PlCn-x  +  P2C2_i) 

<  - 2 -  =  C«- 

[l-a(ac„_1  +  7c2_1)] 

Therefore  the  sequence  {c„}  becomes  a  decreasing  sequence  with  cn  <  1  for  all  n. 
If  0  <  s  <  1  and  cn  <  sc„_i,  then 

_  Pycl  ( P0  +  P\Cn  +  P2c2) 

"+1  [1  -  a(acn  +  7C2)]2 

<  s2Pycl_  1  (Pq  +  PlSC»-i  +  P2s2c2_!) 

[l  -  a  (asc„_i  +  7S2c2_!)] 2 
2  Pvcn-1  (Po  +  PlCn-1  +  P2C2_ 1)  2 

<  S  - - - 2 - “  =  S  Cn- 

[l  —  a  (ac„_i  +  7cn-i)] 

Let  (  =  —  =  Ci  =  ai6i,  then  we  have  0  <  (  <  1  and  Ci  <  £co  =  (,  so  that 
Co 


Cl 

< 

Cco, 

C2 

< 

C2ci, 

C3 

< 

C22c2, 

C4 

< 

>23 

C  c3, 

Cn+1 

< 

£(2n+2 
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On  the  other  hand,  with  the  sequence  {dn}  under  the  assumption  that  aibi  <  1 
we  have 

dn  =  =  ( aCn + 'ycl)  ^ 

<  ( acn  +  7 cl)  ±=acn  +  7 c2 

<  (a  +  7  )c„ 

<  (a  +  7)C2"  •  l 

since  {an}  is  an  increasing  sequence,  and  ao  >  1. 

We  have  thus  proved  the  following  estimates. 

Lemma  3.  We  assume  that  aq&i  <  1.  Then  the  sequence  {cn}  is  a  decreasing 
sequence  and  for  all  n  €  N  we  have  the  following  estimates 

cn+i  <  en+1  ■  b 
dn  <  (a  +  7)C2"  •  ^ 

where  0  <  £  =  a\bi  <  1. 


Lemma  4.  The  sequence  dn  >  0  is  a  convergent  sequence  and  its  limit  is  0. 


Proof.  Since  an  >  1  is  increasing,  then  l/an  <  1  is  a  decreasing  sequence  and 
further  0  <  l/a„  <  1.  Therefore,  an  is  convergent  (it  is  monotonic  in  a  compact 


set)  to  a  limit  L.  Since  dn  =  —  ( - 

&  \Q"n  &n+ 1 

a(L-L)  =  0. 


,  then  dn  is  convergent  to  the  limit 

□ 


Remark  1.  A  similar  approach  would  work  for  some  of  the  lemmas  in  the  paper  [7], 
as  well.  Some  of  their  results,  like  Lemmas  4—7  can  he  simplified  using  a  similar 

OO 

approach:  for  instance,  in  Lemma  7  of  [7] ,  it  is  claimed  that  di  <  00,  but  that 


is  immediate,  since  di  =  lim  di  =  lim  —  (  1  — 

^ J  n,—>oo  C- — ^  n—too  Cl  ' 


i=0 


i= 0 


i—0 

=  -(1  —  L),  where  L 
a 


would  be  the  (finite)  limit  ofl/an. 


Now,  we  state  the  semilocal  convergence  of  the  method  defined  by  (2.2). 

Theorem  5.  Let  X,Y  be  Banach  spaces  and  F  be  a  twice  Frechet  differentiable  in 
an  open  convex  domain  PL  of  a  Banach  space  X  and  BL(Y,  X)  be  the  set  of  bounded 
linear  operators  from  Y  into  X .  Let  us  assume  that  To  =  F'(xo)-1  6  BL(Y ,  X) 
exists  at  some  Xq  €  and  the  following  conditions  hold: 

(!)  \\F'(x)-F'(y)\\<k1\\x-y\\,  x,yePL 

(2)  ||F"(z)||<M,  xen, 

(3)  ||r0||  <  p, 

(4)  jjroF(xo)||  <  7. 

Let  us  denote  a  =  kif3r).  Suppose  that  £0  is  chosen  so  as  to  satisfy  a(a  +  7)  <  1 
and  a\bi  <  1,  where  a  and  7  are  given  by  (2.3).  Then,  if  B(xo,r if)  C  PL,  where 
r  =  Y^-o  dn,  then  the  sequence  {xn}  defined  by  (2.2)  and  starting  at  Xq  converges 
to  a  solution  x*  of  the  equation  F( x)  =  0  .  In  this  case,  the  solution  x*  and 
the  iterates  xn  belong  to  B{xo,rrf),  and  x*  is  the  only  solution  of  F(x)  =  0  in 
B(xq,  2/ (fci/3)  -  rrf)  D  PL. 


Furthermore ,  the  error  bound  on  x*  depends  on  the  sequence  {dn}  given  by 

OO  /  \  oo 

\\xn+i  -  a:*  II  <  dky  <  7  V  ^  C2*,  C  =  aih-  (3.1) 

k—n-\- 1  ^  k=n+ 1 

Proof.  It  is  easy  to  see  that  the  sequence  {xn}  is  convergent.  Hence,  there  exists  a 
limit  x*  such  that  linin^oo  xn  —  x* .  The  sequence  {an}  is  bounded  above  since 


1  —  a  ^2  d  i  1  —  a  '^2  d  i 

2=0  *=o 

Since  limn^oo  dn  =  0,  by  (2.7),  we  have  linin^oo  bn  =  0.  This  indicates  that 
liuin^oo  Cn  =  0.  Thus,  by  (2.26)  and  by  the  continuity  of  F,  we  proved  that 

\\F(x*)\\=0. 


\\xn+l  -  Soil  <  ||®n+l  ~  Xn\\  +  \\xn  —  Xn_i||  4 - h  11^1  ~  ^0 II 

n 

<  ^2  <  r??,  (3.2) 

k- 0 

where  r  =  'f2^=0dn.  We  conclude  that  xn  lies  in  B(xo,rrf)  and  taking  limit  as 
n  — >  oo  we  have  x*  €  B(x o,  rrf). 

To  show  the  uniqueness  of  the  solution,  let  y*  €  B(x o,  2/(fci/3)  —  rrf)  be  another 
solution  of  F{x)  =  0.  Then 

0  =  F(y*)-F(x*)=  f  F' (x*  +  tfy*  —  x*))dt(y*  —  x*).  (3.3) 

Jo 

To  show  that  y*  =  x* ,  we  have  to  show  that  the  operator  fQ  F'(x*  +  t(y*  —  x*))dt 
is  invertible.  Now,  for 

llToll  f  \\F\x*+t{y*-x*))-F\x0)\\dt 
Jo 

<  kiP  f  \\x*  +  t(y*  —  x*)  —  xo\\dt 

Jo 

<  kifJ  [  ((1  -  t)\\x*  -  x0||  +  t\\y*  -  x0|| )dt 

Jo 


-  ~2  y11  +  kf/3  ~  rV)  =  (3'4) 

it  follows  from  Banach’s  Theorem  [1]  that  the  operator  f*  F'(x*  +t(y*  —  x*))dt  has 
an  inverse.  Therefore,  y*  =  x* . 

For  every  m  >  n  +  1,  we  have 

||xm-xn+i||  <  ||xm  -  xm_i||  +  ||xm_i  -  xm_2||  H - b  ||xn+2  -  xn+i  || 

m—  1 

<  ^2  dkr)  <  rrj.  (3.5) 
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By  taking  in  — >  oo,  we  get 


OO 

K+1-*1<  E  dkr]  <  rrj,  (3.6) 

k=n+ 1 

and  from  Lemma  3 

OO  /  \  OO 

\\xn+1  -  z*||  <  E  dkV  <  a  >7  V  E  C2‘,  0<C<1,  (3.7) 

k=n+ 1  k—n-\- 1 

which  shows  that  {:rn}  converges  and  proves  (3.12).  This  completes  the  proof.  □ 


4.  Examples 

In  this  section,  we  give  some  examples  to  illustrate  the  previous  convergence 
result. 

Example  4.1.  [7]  Let  X  =  (7(0,1]  be  the  space  of  continuous  functions  on  [0,1] 
and  consider  the  integral  equation  F(x )  =  0,  where 


F(x)(s)  = —l  +  x(s)  +  Ax(s)  f 

Jo 


-x(t)dt, 


(4.1) 


/  o  s  +  t 

where  s  €  [0,  l\,x  €  (7[0, 1]  and  0  <  A  <  2.  The  norm  is  taken  as  sup-norm.  We 
easily  find 

r1  s  r1  s 

F' (x)u(s)  =  u(s)  +  Au(s)  /  - x(t)dt,  +  Ax(s)  /  - u(t)dt,  u  £  fl,  (4.2) 

Jo  s  +  t  Jo  s  +  t 


and 


F"(x)(uv)(s)  =  Au(s)  /  - v(t)dt  +  Av(s) 

Jo  s  +  t 


s  +  t 


i(t)dt,  u,v  £  Q.  (4.3) 


Since 

III-F'W  -  -F'(2/)](w)||  <  2Aln2||«||||x  -  y||, 
we  get  k\  =  2Aln2.  Since 

f*i  „ 


(4.4) 


|E"(a;)||  <  2Amax0<s<i 


s  +  t 


-dt 


we  get  M  =  2Aln2. 

We  also  have  from  [7]  that  /?  = 


and  rj  = 


=  2A  In  2,  (4.5) 

Iko  -  1||  +  Aln2||a;o||2 


In  our  case,  we  have  a  =  k\f3r)  = 


l-2Aln2||x0||  '  l-2Aln2||x0|| 

2Aln2(||a’o  -  1||  +  A  In  2||cco||2) 


(1  —  2A  In  2||xo||)2 


\q2  ^  _|_  M  _  Q\ 

If  Xq  is  chosen  so  that  a(a  +  7)  <  1,  where  a  =  - - ,7  = 

oz 

A  In  2( ||a;o  —  111  +  A  In  2|| a?o II 2 )  r  .  .  . 

- - — — — — - — - ,  the  sequence  {xn\  dehned  by  (2.2)  and  starting  at 

(1  -  2Aln2||;ro||T 

such  an  xq,  converges  to  a  solution  x*  of  the  equation  F(x)  =  0. 
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Example  4.2.  Consider  the  solution  of  the  nonlinear  equation  F(x)  =  x3  +x  —  10 
on  [1,3].  We  let  9  =  —2.  Now  the  initial  condition  is  Xq  =  1.7,  and  it  is  easy 
to  show  that  F'[x o)  =  3x(j  +  1  =  9.67,  a  =  1,  /?  =  ||E'(a;o)_1||  =  .1034126163, 
r]  =  ||F,(xo)-1F(x0)||  =  | (.  1034126163) (—3.387) |  =  .3502585314  and  ||F"(x)||  = 
||6x||  <  18  =  M.  Now  F'(x)  —  F'(y)  =  3(x  —  y)(x+y)  and  therefore  k\  =  18.  There¬ 
fore  a  =  kifhj  =  .6519807199  and  7  =  f/3y  =  9(.1034126163)(.3502585314)  = 
.3259903600.  Then  the  condition  holds:  a(a  +  7)  =  0.8645201495  <1.  As  a  result, 
the  convergence  of  this  nonlinear  equation  can  he  studied  by  Theorem  5. 


Remark  2.  One  can  take  a  larger  interval,  i.e.  [1  ,A]  for  A  >  3  and  still  satisfy 
the  condition.  Suppose,  we  let  x$  =  A  >  3,  then  F'(A)  =  3 A2  +  1,  a  =  1,  (3  = 

ll^,(^o)_1||  =  pp+r ,  V  =  ll^'(£Co)_1^(a;o)||  =  AIa2+i°  and  \\F"(X)\\  =  IIM  < 
6^4  =  M.  Now  F'{x)  —  F'(y)  =  3(x  —  y)(x  +  y)  and  therefore  £7  =  6A.  Therefore 
a  =  kiflr)  =  6A^3J~2^_~|P  and  7  =  =  SA^+f^® .  Then  the  condition  holds: 


n((y  sy)  —  fi  4 -A  +A—10  1  o  A  A  +A— 10 
a\a  -1-  1)  —  O/i  (3^2+1)2  1  -h  o/i  (3^2+1)2 


t  ^  1  as  A 


00. 


Clearly  the  number  of  iterations  (say,  n)  required  for  convergence  depends  on 
how  close  A  to  x*  =  2.  Experimenting  with  various  values  of  xg  yields  the  following 
results: 


x0 

=  1.7 

n 

X0 

=  3. 

n 

X0 

=  5. 

n 

X0 

=  10. 

n 

Example  4.3.  Let  us  consider  the  system  of  three  nonlinear  equations  F(x,  y,  z)  = 
0  where  F  :  LI  — ►  1Z 3  where  14  =  [07,61]  x  [02,62]  x  [03,63]  is  a  domain  of  F 
containing  a  solution  of  this  system,  and 

F(x,y,z)  =  (x2+y2+z2—4,,x2+y2+z2  —  2x—2y—2,x2+y2+z2—4:y—2z—A).  (4.6) 


Then  we  have 


F'{x,y,z ) 


2x  2 y  2 z 

2x  —  2  2r/  —  2  2  z 

2x  2y  —  4  2z  —  2 


(4.7) 


and 


F'(x,y,  z)-1 


1 

2(— y  +  2  z  +  x) 


—y  +  z  +  1  y  —  2z  z 
z  —  1  +  X  —x  —z 
— x  —  y  +  2  2x  —x  +  y 


(4.8) 


We  use  the  Frobenius  norm  inTZ?:  ||X||  =  (x2  +  y2  +  22)1/2  for  X  =  ( x,y,z )  €  TZ3. 

/  3  3  4/2 

The  corresponding  norm  on  A  £  1Z3  x  7 Z3  «PII=  E-=iEi=iKf  • 
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The  second  derivative  is  a  bilinear  operator  on  7Z3  given  by 


flxx 

flxy 

flxz 

flyx 

foyy 

flyz 

flzx 

fozy 

fizz 

f2xx 

foxy 

f2xz 

foyx 

foyy 

foyz 

J2ZX 

fozy 

j2zZ 

f?>xx 

fsxy 

f?>xz 

foyx 

foyy 

fsyz 

fszx 

fozy 

f?>zz 

2 

0 

0 

0 

2 

0 

0 

0 

2 

2 

0 

0 

0 

2 

0 

0 

0 

2 

2 

0 

0 

0 

2 

0 

0 

0 

2 

For  a  bilinear  operator  B  :  1Z3  x  1Z3  1Z3 ,  given  two  vectors  X,  Y  €  7 Z3,  we  have 


'  B1  ' 

2/i 

B(X,Y)  :=  (xi,X2,x3) 

b2 

2/2 

.  Il‘‘  . 

_  2/3 . 

b\lx  i  +  b\2X2  +  bfx  3 
bfx  i  +  bfx2  +  bfx3 
bfx  i  +  bfx  2  +  bfx  3 


bfx  i  +  bfx2  +  bf  x3 
bfx\  +  bfx2  +  bfx3 
bfxi  +  bfx  2  +  bfx3 


bfxi  +  bfx2  +  bfx3 
bfx\  +  bfx  2  +  bfx3 
bfx  i  +  bfx  2  +  bfx3 


where 


fe}1 

bf 

bf 

bf 

bf 

bf 

bf 

bf 

bf 

Bx  ' 

— 

- 

— 

bf 

bf 

bf 

= 

bf 

bf 

bf 

— 

bf 

bf 

bf 

S3 

— 

- 

bf 

bf 

bf 

bf 

bf 

bf 

bf 

bf 

bf 

(4.10) 


(4.9) 


y  i 

2/2  , 
2/3 
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We  consider  the  norm  of  a  bilinear  operator  B  on  TZ?  by 


\\B\\ 


sup 

IHNi'V 


EEEft 

i= 1  j=l  \fe= 1 


Jk 


Uk 


(4.11) 


where  u  =  (ui,v,2, 113). 

In  our  case,  for  any  triple  ( x,y,z ), 


M=\\F"(x,y,z)\\ 


Now,  since 


2V3  sup  {|mi|2  +  M2  +  |U3|2}1/2  =  2V3.  (4.12) 
IMI=i 


x  — 

u 

y 

—  V 

z  —  W 

F'(x, 

y,z )  - 

F'(u,v , 

w)  =  2 

x  — 

u 

y 

—  V 

z  —  w 

x  — 

u 

y 

—  V 

z  —  w 

we  have 

II  F'(x 

,y,z) 

-F\u 

> v,  u?)  || 

=  2a/3  [(x 

- 

uf 

+  (y 

-v)2  + 

II 

to 

y> 

z)~ 

-  (u, 

v,w)  II, 

(4.13) 


w) 


21 1/2 


and  therefore  k\  =  2y/3. 

We  let  8  =  —2,  and  so,  a  =  1.  If  we  choose  Xo  =  —  0.2,  yo  = 
1.2,2o  =  -2.1,  then  /3  =  \\F'(x0,  y0, 20)_1||  =  .675981901098132166,??  = 
\\F' (xo,  yo,  Zo)~1F(xo,yo,  Zq)\\  =  .28125.  Therefore  a  =  ki(3r]  =  .6585946863  and 
7  =  4t/3?7  =  .3292973430.  Then  the  left,  side  of  the  condition  holds  a(a  +  7)  = 
.8754681666  <1.  As  o  result,  the  convergence  of  this  system  of  equations  can  be 
studied  by  Theorem  5. 
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